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We investigate the effect of chemical freeze-out on identified particle spectra at 200AGeV 
Au-Au Collisions at RHIC, by utilizing a full three-dimensional hydrodynamical calcula- 
tion. The hydrodynamical code SPheRIO we employed is based on the smoothed particle 
hydrodynamic algorithm. In order to describe the spectra of strange hadrons, the code 
has been further improved by explicitly incorporating the strangeness conservation and 
a chemical freeze-out mechanism. In our model, strange hadrons such as A, S, Q and <j> 
undergo the chemical freeze-out immediately after the hadronization, and their multi- 
plicities are fixed thereafter. At a lower temperature the thermal freeze-out takes place 
for all the particles. It is shown that the present model provides a reasonably good de- 
scription for the spectra of identified particles, in particular, considerable improvement 
is observed for those of strange hadrons. 



1. Smoothed particle hydrodynamic model 

Statistical model and hydro-inspired model calculations indicate the picture of 
early chemical freeze-out(CFO) for strange hadrons in relativistic heavy ion 
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collision !™™ The CFO temperature can be extracted by fitting the model 
calculation of particle ratios for hadrons to the experimental data. On the other 
hand, one obtains the thermal freeze-out temperature from the slope of transverse 
momentum distributions by assuming some radial flow profile. Naturally, it is there- 
fore motivating to see how this scenario can be adopted by more realistic models, 
where a full three-dimensional hydrodynamical evolution is taken into account. In 
this work, we present results of hydrodynamical model calculations, inspired by this 
picture. The hydrodynamica l mo del we employed is based on smoothed particle hy- 
drodynamic(SPH) algorithm-'- 12 ', In this model, the matter flow is parametrized in 
terms of discrete Lagrangian coordinates, of the so-called SPH particles. As a result, 
the hydrodynamic equations are reduced to a system of coupled ordinary differen- 
tial equations. The code which implements the entropy representation of the SPH 
model for relativistic high energy collisions, and which has been developed within 
the Sao Paulo - Rio de Janeiro Collaboration, is called SPhcRIO. As it has been 
shown, the model is an efficient and robust method to tackle the problems concern- 
ing relativistic high-energy nucleus-nucleus collisions, which are characterized by 
highly asymmetrical configurations. It has been successfully used to investigate the 
effects of the initial-condition fluctuations and a dopting the c ontinuous emission 
scenario for the description of decoupling process | 8 | 7 | 9 | 10 | 11 | 12 | 

In the present calculation, the model has been further improved in order to con- 
sider strangeness conservation and to adopt the scenario of CFO. The strangeness 
conservation is implemented by explicitly incorporating strangeness chemical po- 
tential into the code, and correspondingly a different set of equation of state(EOS) 
has been built and utilized. As a good approximation, we assume local strangeness 
neutrality throughout the hydro evolution. 

To adopt the picture of CFO, as a first step, we take a simple postulate that 
for strange hadrons such as A, S, f2 and <fc, the CFO takes place immediately after 
they are completely hadronized. These particles cease to have inelastic collisions 
and therefore their abundances are fixed and determined only by partonic EOS. 
Instead of introducing a specific temperature, CFO is incorporated in a simple and 
parameter-free approach. This is inspired by the earlier estimation from statistical 
mode0that CFO takes place close to the phase transition region. It is also because 
at this point, we would like first to check whether the picture of CFO would quali- 
tatively improve the result of hydrodynamic model without going into the details of 
parameters. With the system being cooled and rarefied further, the thermal freeze- 
out occurs at a lower temperature Tf which serves as an adjustable parameter in 
our calculation. 

2. Results and discussions 

As in the previous calculations, we use NEXUS event generator to produce event- 
by-event fluctuating initial conditions. We make use of a rescaling factor-^ to fix the 
pseudo-rapidity distribution for all charged particles. In the present calculation, we 
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introduce another parameter a which provides initial transverse velocity in addition 
to the one provided by NEXUS. The reason for the introduction of a is because 
NEXUS almost does not provide any initial transverse velocity. Thus the initial 
transverse expansion reads 



v T = NEXUS + tanh(or)f 
where r is the radial distance from orgin. 



(1) 



exp from PHOBO 
- EOS with strangeness 




Fig. 1. The pseudo-rapidity distribution for all charged particles for the most central Au+Au 
collision at 200A GeV. The data are from PHOBOS Collab. 



As a preliminary application, we have calculated the pseudo-rapidity distribu- 
tion for all charged particles. In Fig.l, we present the pseudo-rapidity distribution 
for the most central Au+Au collisions at 200A GeV. The experimental data are 
from PHOBOS Collab, taken in the most central Au+Au at 200A GeV. 

We show in Fig. 2, the experimental transverse momentum distribution data for 
all charged particles, which can be well reproduced with a choice of a = 0.04fm _1 
and freeze-out temperature Tf = 110 MeV. The experimental data are from STAR 
CollarJSl, taken in the most central Au+Au at 200A GeV, with r) = 0. 

With the parameters chosen, in the following, we calculate the spectra for var- 
ious hadrons. In Fig. 3 we show the transverse mass spectra of pions, protons and 
kaons for most c entr al collisions at mid rapidity, as well as experimental data from 
BRAHMS CollatM The same spectra for A, S and f2 are depicted in Fig. 4, to- 
gether with data from STAR CollalM We use dotted lines to represent the results 
obtained without incorporating CFO, and solid lines for those with CFO switched 



February 9, 2008 0:57 WSPC/INSTRUCTION FILE 

4 Wei-Liang Qian 



wlqian 



10 4 i 



IO- 
CS 



A STAR Data, ti=0 

Event by Event, T=1 1 0MeV 



T3 

qT 10"' 



Vf-i 




10'*- 



P T (GeV) 



Fig. 2. Transverse momentum distribution of all charged particles for the most central Au+Au 
collision at 200A GeV in the rapidity interval —1.0 < r\ < 1.0. The data are from STAR Collab. 




P T (GeV) 



Fig. 3. Transverse momentum distribution of pions, kaons and protons for the most central Au+Au 
collision at 200A GeV in the rapidity interval —1.0 < r\ < 1.0 for the most central collisions. The 
data are from BRAHMS Collab. 



It is observed that the present hydrodynamic model gives good description of 
the experimental transverse momentum spectra for pions, kaons and protons even 
without turning on CFO for these particles. Although it gives the correct slopes for 
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Fig. 4. Transverse momentum distribution of A, H and f! for the most central Au+Au collision at 
200A GeV in the rapidity interval -1.0 < n < 1.0. The data are from STAR Collab. The dotted 
lines indicate the results obtained without CFO, and solid line for those with CFO incorporated 
in the model. 

the spectra of strange hyperons such as A, 3 and f2, the disagreement comes from 
the multiplicities of the spectra. By introducing CFO, those results are improved 
significantly. While the spectra for pions, kaons and protons almost remain the 
same, (this is not shown in the figures.) it provides a good fit of strange hyperons 
A, H and 0. The CFO amplifies the multiplicities of strange hadrons with respect to 
the ones in the original model owing to higher CFO temperature, as compared with 
the thermal freeze out one. Meanwhil e th e slopes of the spectra remain unchanged. 
It was indicated experimentally in RefP, as one goes to large rapidity region where 
the baryon density differs sizably from zero, the outcome of CFO would be much 
more significant. Further work on this topic is under progress. 

It is worth no ting that the effect of CFO has been discussed by several authors. 
In reflHIinEQEI] the early CFO is studied in term of hybrid model where the hy- 
dodynamic evolution is complemented with a hadronic cascade model. The pure 
hydrodynamic evolution calculation has been carried out either only in the trans- 
verse plane ass umi ng Bjorken's scaling^ 1 or with the baryon chemical potential 
taken to be zerc^l. 

We acknowledge financial support by FAPESP (2004/10619-9, 2005/54595- 
9, 2004/13309-0 and 2004/15560-2), CAPES/PrOBRAL, CNPq, FAPERJ and 
PRONEX. 
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